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1 Introduction 



\ The metric problem of equivalence is the problem to decide whether two pseudo- 

' Riemannian spaces are locally isometric. This question is of importance in general 

relativity ([39, Ch. 9]), where isometric means gravitationally identical. 
\ Existing classifying algorithms depend on computation of curvature invariants, 

0^ ' which are components of Riemann curvature tensor and its covariant derivatives 

' in a particular frame specified in the course of the algorithm. This is the central 

. idea of Cartan's [8, 6] solution to the equivalence problem by the method of moving 

frames, and the Karlhede [22, 24] algorithm derived from it. Karlhede classification 
lies behind two projects of publicly available databases of exact solutions of Ein- 
stein equations [21, 36]. In vacuum case, the Riemann tensor can be replaced with 
\ its traceless part, the Weyl tensor. Presently, the Cartan-Karlhede curvature invari- 

ants are the only invariants known to solve the metric equivalence problem in full 
generality. 

The main difficulty connected with the Cartan-Karlhede algorithm lies with the 
size of the data generated at the nth iteration of the main loop. The number and 
complexity of generated curvature invariants grow quite rapidly with the order n of 
covariant derivatives. During the past two decades, much effort has been exerted to 
tight the bounds on n, with recent limit values being collected in [39, Ch. 9, Ta- 
ble 9.1]. The required n can be as high as five, even though the upper limits are 
reached only in very special cases [38]. Understandably, the higher n is, the bigger 
portion of generated invariants are actually redundant even if algebraic identities are 
accounted for as in MacCallum and Aman [30]. This is the price paid for leaving 
no metrics aside. From the computational perspective, one is naturally interested in 
solving the equivalence problem in terms of minimal possible number of invariants of 
lowest possible order. In this paper we provide such a solution within the class T of 
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generic metrics that possess two commuting orthogonally transitive Killing vectors. 
Our choice reflects the abundancy of exact solutions of class T in general relativ- 
ity, be it stationary axisymmetric [10] or cylindrically symmetric spacetimes. The 
orthogonal transitivity is an additional simplification that not only occurs for all 
explicitly known stationary axisymmetric metrics [39, p. 294], but can be derived 
from much simpler property of invertibility [35, 9]. A fairly large portion of explicitly 
known spacetime metrics belong to this class, mainly due to the fact that powerful 
generation techniques operate within T or its vacuum subclass. Reduction to a rather 
compact Ernst cqiiation [17] was the key to solution by the inverse scattering method 
of Belinsky and Zakharov [4], as well as a number of Backlund transformations, such 
as the Harrison transformation [19], Hoenselaers-Kinnersley-Xanthopoulos trans- 
formation [20], Kinnersley-Chitre transformation [28], and Neugebauer transforma- 
tion [31]. The majority of these results have been extended to metrics with suitable 
energy-momentum tensors, notably electro-magnetism. For relations between differ- 
ent solution-generating techniques see Cosgrove [14]. For a number of explicit solu- 
tions see [39, 40]. As shown by Kinnersley [27], every vacuum Petrov type D metric 
belongs to T, the result having been extended to aligned electro- vacuum metrics with 
cosmological constant by Debever and McLenaghan [15]. On the other hand, metrics 
from the class T can be of any Petrov type except ///. Last but not least, Cos- 
grove [12] almost solved the equivalence problem for vacuum metrics in this class. 
It should be noted that metrics of class % pass the Karlhede classification rather 
smoothly, requiring no more than two covariant derivatives to be taken. The Rie- 
mann tensor itself being of order two, the computed invariants are of order up to four 
in terms of metric coefBcients. It may therefore come as a surprise that appropriately 
chosen invariants of first and second order suffice, at least for generic metrics. Our 
goal is to demonstrate that. 

The content of the paper is as follows. In Section 2 we review the transformations 
that preserve the class T. Then we compute the number of lowest order invariants 
to show that four of them are of the first order. 

In Section 3 we present the four independent invariants Cp,Cy., Q^, as well as 
an invariant frame X, y on . They provide us with a solution to the equivalence 
problem rather immediately. It should be stressed that we consider only generic met- 
rics within our class. The exact meaning of "generic" is specified along the solution of 
the equivalence problem. In particular, metrics are required to satisfy the condition 
that their Killing algebra possesses a unique two-dimensional commutative subalge- 
bra. However, we leave open the problem of recognizing this condition in terms of 
first-order invariants. 

Our work was to a great extent inspired by Cosgrove's treatment [12] of the 
subclass of vacuum space-times. Exploiting an auxiliary invariantly defined metric of 
constant negative curvature, Cosgrove [12, 13] devised a method to find new vacuum 
solutions of Einstein equations without cosmological constant. The work [12] solves 
the equivalence problem for empty spacetimes modulo transformations of the natural 
harmonic coordinate r and its companion coordinate. 

In Section 4 we attempt a comparison between our and Petrov/Karlhede classi- 
fication. The Cosgrove invariant is rather immediately identified with imaginary 
part of ^'2, another, C^, with the real part of ^'2 in case of Ricci-flat (i.e., vac- 
uum) metrics. We were unable to identify the others. The last section is devoted to 
examples and discussion. 

Beyond the scope of this paper we leave investigation of the non-generic cases as 
well as extension to wider classes of metrics. A general classification of conditions that 
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can be put on a two-dimensional Killing distribution and its orthogonal complement 
can be found in [11, 3] along with the description of general solution in nongeneric 
cases. 

2 The metric and the pseudogroup 
2.1 The metric 

In this section we introduce the class T of pseudo-Ricmannian metrics of interest in 
this paper and the pseudogroup of local diffcomorphisms they admit. Let be a 
four-dimensional manifold. We consider a non-degenerate metric g on possessing 
a unique two-dimensional commutative algebra ^ of Killing vectors. Let S denote the 
vector distribution on A4 generated by ^. Then H is two-dimensional, since no non- 
constant multiple of a non-zero Killing vector can be a Killing vector. Let ^i, ^2 denote 
arbitrary generators of the algebra A. We assume additionally that the Killing vectors 
be orthogonally transitive [35, 9], meaning that the distribution S"*- generated by the 
covectors g(^i, ■ ) is Probenius integrable and transversal to H. Hence the tangent 
spaces decompose as TaAi = © "^a- By Frobenius integrability, the manifold 
becomes locally foliated by a two-parameter family S of two-dimensional leaves of 
the distribution S-*-. For an arbitrarily chosen leaf Sq, points of <So parametrize the 
orbit space at least locally. Note that global properties of the foliation are 

irrelevant to the local equivalence problem. 

Metrics of class T have a simple description in terms of adapted coordinates (see 
also Weyl-Lewis-Papapetrou coordinates below). Denoting by »So a fixed surface of 
the family 5, wc arbitrarily choose coordinates t^,t'^ in an open subset of Sq and 
extend them along the trajectories of the fields to functions in an open 

subset U of M.. Shrinking U if necessary, we determine functions u^,u^ in U by the 
requirement that ^lu'^ = and SoDU = {u^ = = 0}. In the coordinate system 
{t^,t'^,u^,u^) we have = d/du'^, k = 1,2, while the metric g assumes the form 

g = gij df dt^ + hki du'' du\ (1) 

where gij = gji, hki = hik are functions on S, i.e., depend on t^,t'^ only. The induced 
metric gj^ on S \s g = gij df dP . Obviously, 

det g = det g det h ^0. 

The other leaves of the foliation S are = consti, = const2. Since shifts along 

the trajectories of Killing vectors are isometrics, the particular choice of the leaf Sq 
is irrelevant. Henceforth we drop the index and write simply S instead of Sq. 

The metric g can be identified with a section of the second symmetric power 
S'^T*M. — > Moi the cotangent bundle T*M — > M. The two-dimensional spaces S^, 
a G 5, constitute a vector bundle H — > S. Considering the decompositions TaM = 
'EaXTaS for a £ S, the metric g = g\s can be identified with a section of S'^T*S — > S, 
while the remaining metric coefficients h^i = g(^fc,^;) parametrize sections of the 
bundle S^E* — > S. Thus, locally, metrics g are in a one-to-one correspondence with 
sections {g, h) of the six-dimensional vector bundle vr : S'^T*S x S'^E* — > S. This 
correspondence, given simply by (1) in adapted coordinates, is only valid under the 
condition that the distributions H and H"*- are fixed. 
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2.2 The pseudogroup 

Turning back to our classification problem, we need to know the corresponding local 

transformation pseudogroup acting on sections of the bundle tt. Otherwise said, what 
are local diffeomorphisms $ of such that g' = <I>*g is of the same symmetry 
properties as g (both have the same algebra hence the same generators ^i, ^2 in 
and induce one and the same decomposition TaM = S^)? More to the point, 
what they are in terms of g and /t? 

Let <3 denote the pseudogroup of local diffeomorphisms of the base S of vr. These 
have a unique lift along the bundle S^T*S — > S, given by g' = (j)*g, 4> £ &. Next 
consider the standard representation of GL2 in the two-dimensional vector space of 
the algebra ^. That is, given (af) G GL2, let .^^ = a^,^.j. Then we have the induced 
representation in each H^, hence in each S'^H*, namely /i^^ = = ^k^j^^ij- 

Thus we arrive at the product <5 x GL2, henceforth referred to as the pseudogroup 
of transformations of the total space S'^E* x S'^T*S of tt. 

Proposition 1 Let two metrics g and g' on M. possess a unique two-dimensional 
commutative algebra R of Killing vectors, which induces one and the same decomposi- 
tion TaM = Ha©H;|-, a £ M. Then g' = $*g locally with respect to a dijfeomorphism 
^ of Ai if and only if the corresponding sections (g, h) and (g', h') of ir are locally 
equivalent with respect to the action of & x GL2 given above. 

Proof. To prove the 'if part, let (j) denote a local diffeomorphism on S such that 
g' = (p*g and let (a*) G GL2 be such that /i^^ = a\ajhij. Denoting by $ the local 
diffeomorphism of A4 with coordinate description = a^u\ o # = o ^, from 

(1) we get $*g' = g. 

To prove the converse, consider a local diffeomorphism $ of such that g' = $*g. 
Let ^1, ^2 be the basis of ^. Then <5*^2 arc Killing vectors for g and belong to 

the unique commutative subalgebra ^, therefore they are linear combinations of ^1 , ^2 
with constant coefficients. Denoting by a the local diffeomorphism of M. with local 
description o a = afu^ (a^.) G GL2, f o a = f, we have a^^fc = Hence, 
$ o preserves both vectors ^i,'^2- The coefficients (a^ G GL2 being uniquely 
determined, a can be regarded as a representation of GL2 by local diffeomorphisms 
of M, equivalent to the standard representation of GL2 in each S^. 

Now, ^a~^S = #5 is another leaf of the family, i.e., I3S, where /3 is a shift along 
the trajectories of ^i,^2j i-e., o /? = + Denoting (f> := (3~^ 0^0 a~^, we 
sec that is a local diffeomorphism that preserves S as well as both vector fields 
^1,^2- Moreover, cf) is uniquely determined by its restriction hence effectively 
belongs to the pseudogroup (5. Thus the local diffeomorphism $ can be decomposed 
as^ = (3o(j)oa into an automorphism a G GL2, a local diffeomorphism of (j) £ 
and a shift. However, shifts are isometrics, hence g' = $*g = a*4>*P*g = a*(f)*g. 
Thus we are left with the pseudogroup (5 x GL2. □ 

The classification problem for metrics g thus reduces to identifying orbits of the 
action of x GL2 on sections of the bundle tt. Following Lie's classical method, 
one has to find a sufficient number of independent scalar differential invariants of 
the action. These are functions on the jet prolongation J°°tt that are invariant with 
respect to the infinitesimal action of © x GL2 on J°^tt. Thus the first steps to be 
done are to describe the action of the corresponding Lie algebra g x 0I2 on tt and its 
extension to the jet prolongation J°°7r. 
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Starting from the description of elements of © as diffeomorphisms (f) given by 
o (j) = (f)^ {i} ,t^) ^ we see that q acts by vector fields 

^^ = V''|l (2) 

with ijj = {tjj^^t^ ,t'^),'tjj'^{t^ ,t'^)) being an arbitrary couple of functions on S. To 
determine the extension of infinitesimal transformations (2) on tt, we require 
that the Lie derivative U^{g) = U^{gij df dP) be zero, i.e., 

__dr_ or 

'^tpdlj — 9lS r^j 9sj ■ 

Thus we have 



The remaining infinitesimal generators are Va with a G corresponding to 
action of GL2 on components hki- In terms of coordinates u'' on M we have Va = 
aliu\d/du^) with a = {a^) G being an arbitrary constant 2x2 matrix, hence 

Va^ki = —hksdi — hsial 
similarly as above. Then 

K = - J2(^ksat + hsial)^, (3) 
k<l 

since components transversal to S are suppressed on tt. 

Next we consider the infinite-dimensional vector bundle J°°tt. Recall that obvious 
coordinates along the fibres of J°°tt are the formal derivatives of Qij and hki of all 
orders with respect to e.g., gij^k = dgij/dt^, etc. On J°°7r one has the usual 

total derivative 

TOO " — TOO 

By := U,^ and Va ■= Va ^ we denote the well-known prolonged fields on 

J°°7r, characterized by the defining relation [U^,Dk] = —{dijj^/dt^)Ds. For details 
and explicit formulas see [1, 5, 32, 33]. As is well known, scalar differential invariants 
can be identified with functions on J°°7r invariant with respect to the fields U^p 

and Va- These functions form a commutative associative M-algebra, which can be 
thought of as algebra of functions on the orbit space J°°7r/(© x GL2). 

We finish this section with a proposition giving the number A^^ of independent 
rth order scalar differential invariants, as found by the classical method of Lie: 

Proposition 2 The dimension Nr of the orbit space J'"7r/(0 x GL2) is given by the 
following table: 



r 





1 


2 


3 


4 ... 


Nr 





4 


14 


28 


46 ••• 
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The proof reduces to routine counting the number of independent equations in 
the system U^f = 0, Vaf = on J^tt. Computations are considerably simpler when 
using the Weyl-Lewis-Papapetrou metric coefficients (see Sect. 3.4 below). 

Remark 1 A comment on a possible source of misunderstanding is due. Proposi- 
tion 2 refers to scalar differential invariants as functions on the jet space J°°7r. If such 
a function, say F, is evaluated for a particular metric g, then it becomes a function 
on S, which we shall denote as F\g (formally _F|g = F o j°°ag, where denotes a 
jet prolongation of a section of the bundle tt and dg is the section associated with g). 
Analogous correspondences hold for other geometric objects such as forms and vec- 
tor fields. Hence another interpretation of scalar differential invariants as functions 
on S. 

Both interpretations arc natural and indispensable. For instance, the order of an 
invariant can only be seen in the context of jet spaces, while the most natural way 
to construct an invariant consists in combining various invariant geometric construc- 
tions on S, as demonstrated in the next section. It is usually harmless to use one and 
the same notation with both interpretations and omit the symbol |g. However, one 
should bear in mind that independence of functions on S is very different from that 
on J°°iT. The maximal number of independent functions is two on S, and unlimited 
on J°°7r. 

3 The equivalence problem 
3.1 First-order invariants 

It is the aim of this section to provide explicit formulas for the four independent scalar 
invariants of the first order predicted in Proposition 2. We utilize invariant geometric 
objects associated with the metric g and with the triple of functions /iii,/ii2 = 
h2i, ^22 on S. Indices are raised and lowered with the metric g. Comma notation is 
used for partial derivatives taken with respect to coordinates t^,t'^ on S. 

Geometric objects on S (functions, metrics) are obviously ©-invariant. There- 
fore, of interest are GL2-invariant geometric objects on S. Of them, GL2-invariant 
functions are the scalar invariants sought. 

The obviously GL2-invariant metric g alone has no invariants of order less than 
two (it has one invariant of order 2, the scalar curvature R, and A; — 1 more indepen- 
dent invariants of order k for each k > 2; see Zorawski [41]). 

However, if a = aij dt^ dP is another GL2-invariant quadratic (meaning bilinear 
symmetric) form on 5, then the trace Cq, = a\ = g^-'aij is a GL2-invariant scalar. 
Moreover, the volume forms vol^ = Vldet g\ dt^ A dt^ and vol^ = \/|detQ;| dt^ A dt"^ 
are a multiple of each other, hence Qa = deta/detg is one more GL2-invariant 
scalar on S. Alternatively, invariants Ca and Qa can be defined as coefficients of the 
characteristic polynomial det(Q; — Xg)/ detg in A. 

Denote by x the function 

X = det h 

on S. Wc already know that x ^ everywhere since otherwise g would be degenerate 
at some point. It is easily checked that x is invariant with respect to the subgroup 
SL2 C GL2, but not to its complement R C GL2 consisting of scalings x 1 — > ex. 
The covector X = dx/x is then invariant with respect to the full GL2 since scalings 
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cancel out. Let p denote the degenerate quadratic form (dx/x)^ = {l/x^)x^iXj df dP 
on S. The trace 



„2 9''^,i^,j 



is a GL2-invariant function on 5, while Qp is obviously zero. 

Another GL2-invariant metric on 5 is x = (dhu d/122 — dhi2 dhi2)/x with com- 
ponents 



Xij 



1 

2x 



h21,i 



hi2,j 

h22,j 



+ 



h21,j h22,i 



Hence two more GL2-invariant functions on S, namely 



X 



hii,i hi2,j 

h21,i h22,j 



det^ 



Next consider linear combinations x+cp, which are GL2-invariant quadratic forms 
on 5 as well. For the particular value c = — | we obtain the Cosgrove form 



1 = X 



This name reflect the fact that 7 coincides with [12, Eq. (2.3)]. While C-y = C~^ — \Cp^ 



the other invariant 

det7 



1 



det g 4x det g 



/ill 


hi2 


h22 




hi2,i 


/i22,l 


hn,2 


hl2,2 


h22,2 



(4) 



is, in general, functionally independent of the invariants Q^, C^, Cp (on J^vr; see 
Remark 1). If so, then we choose Cp, C^, Q^, Q-y to be the basic first-order invariants 
all other first-order scalar invariants are functions of. 

3.2 Invariant frame 

The invariant covector dx/x on S can be turned into an invariant vector by rising its 
index. If done with the help of the metric g, the resulting vector X has components 



X' = g'^ 



; 1 dx 



xdt^' 



If done with the help of the volume form vol^, the resulting vector Y has components 

1 dx 



:V\d^\ dt^ ' 



where e^^ = e^^ = 0, e 



12 



— e 



21 



— 1. It is perhaps worth noticing that X lux = Cp 



while y In X 
Since 



0. 



X' Y' 
X^ y2 



a. 



vTdet 



9\ 



vectors X, Y form an invariant frame at each point of S where Cp ^ 0, i.e., where X 
is not null. Components of g with respect to this frame are 



giX, X) = Cp, giX, Y) = 0, g{Y, Y) = ±Cp, 
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where the sign is that of det g. Components of p are simply 

p(X,X) = C2, p{X,Y) = 0, p{Y,Y) = 0. 

For components of x we easily get 

x{X,X)=CpC^-4Q^ + 4Q^, 
X{Y,Y) = ±4{Q^-Qj), 
x{X,X)x{Y,Y) - x{X,Y)'' = ±ClQ^, 

hence 

X{X,Y) = V±C^pQ^ T <Qx - Qi){CpC^ - 4Q^ + 4Q^) . 
3.3 Second-order scalar invariants 

Lie derivatives of invariant objects with respect to invariant vector fields such as 
X, Y are invariant objects again. Moreover, the vector fields on J°°it corresponding 
to X, Y (see Remark 1) are invariant differentiations, since they commute with the 
vector fields C/^, V a for all V' and a. In particular, if p is a first-order scalar invariant, 
then Xp, Yp are second-order scalar invariants. 

Starting with the four basic first-order invariants Cp,C^,Q^,Q^, we obtain up 
to eight independent second-order invariants in this way. As we shall demonstrate 
below, these suffice for solving the classification problem in generic case. 

Perhaps the best-known scalar second-order invariant of g is the scalar curva- 
ture i?, which is the only independent invariant of second order (in terms of metric 
coefficients) on a two-dimensional manifold like S. 

The invariant vector fields X, Y are orthogonal with respect to g but do not 
commute in general. The coefficients o, b in the commutation relation [X, Y\ = aX + 
bY are second-order scalar differential invariants again. 

Possible ways to construct invariants of second order include finding GL2-invariant 
quadratic forms on 5, such as the Lie derivatives £x X ^^nd £y x, or symmetric 
products dp dq for various choices of first-order scalar invariants p, q. Another 
obvious candidate is the Ricci tensor Ric^j (which, unlike the Ricci tensor Ricjj of 
g, does not necessarily vanish in case of vacuum metrics). Next, associated with any 
smooth function / on 5 is the hessian Hess /, which is a quadratic form defined by 
Uess f{Z,Z')=izVz'df, i.e., 

HesSjj / = f^ij — Tj^jf i^. 

Here V is the Levi-Civita connection associated with g and F^j are the usual Christof- 
fel symbols for V. The trace of Hess / coincides with the Laplace-Beltrami operator 

A/ = g'^ Ressij f. 

Note that the order of Hess / equals two plus the order od /. To obtain a GL2- 
invariant of second order, / must be of order zero. Choosing / to be the logarithm 
of the SL2-invariant x of order zero, we get 

Vij = HesSjj In X = — HesSy x + pij 

(note that scalings of x cancel out). 



8 



Therefore, a number of second-order invariants result as C^, Qa for a running 
through various hnear combinations of Ricjj, fjj, Pij,Xij- (£x x)iji (£y x)ij- 

OnAi, we have of course also the geometric objects associated with the full metric 
g. They can be restricted to S, thereby providing another source of invariant objects. 
Since at most 14 second-order invariants can be independent (as functions on J^vr), 
it should not come as a surprise that, for instance, the restriction of the full Ricci 
tensor Ricjj is a linear combination of the Ricci tensor with f, p and x- 

{'Ric\s)ij = Ricij - \ Vij - \ pij + \ Xij, j = 1, • • 2. (5) 
The remaining components of the full Ricci tensor can be written as 

R'ic2+j,2+j = — 5 ^hij + I Xhij — \ C^hij i, j = 1, . . ., 2. 
The full scalar curvatures are then related by 

R = -R — Cjy — ^Cp-|-^ (6) 
with Cn = Alnx. 



3.4 Weyl— Lewis— Papapetrou parametrization 

The four basic invariants Cp,C-^,Q-^,Q^ look considerably simpler when the coeffi- 
cients hij are expressed via the (three of four) Weyl-Lewis-Papapetrou [29] variables 
r,s,w (r, s nonzero) according to the identification 

h = -[{du^ +wdu^f^s^ {du'^)% 
s 



I.e., 



s 



rw 

s 



rw r{w ^ s 
s s 



(7) 
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in matrix notation. The minus sign corresponds to det h = — < and the plus 
sign to dc t h = r"^ > 0. If is non-null, which can always be achieved, then /in ^ 
and r = \/|det h\ , s = r/hll, w = hi'^/hu. 
Under this parametrization we have 



X = ^r 



^dr^, 



1 



7 



ids'^Tdw'^). 



(8) 



Hence Cosgrove's [12] crucial observation that 7 is a metric of constant curvature 
— 1 if nondegenerate, i.e., if 7^ 0. Formulas 



-9 



1 



=F 



s det 5 



So Wo 



(9) 



will be needed in the sequel. In case of Lorentz metrics we have det g det h <0, hance 

Q-y < 0. 

Proposition 3 // 7^ 0, then the Killing algebra of 7 coincides with the sfe com- 
ponent of the algebra of infinitesimal transformations (3) acting on S. 
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Proof. Both algebras have one and the same set of generators 

dw^ ds dw^ ds 2 dw 
in terms of Weyl-Lewis-Papapetrou variables. □ 

Further simplification is achieved if taking into account that all two-dimensional 
metrics are conformally flat. Then g can be cast in one of the following two forms: 

g = f {d£- dt^ + dt^ df) or g = 2fdt^dt^. (10) 

The full set of Weyl-Lewis-Papapetrou variables contains / as the fourth variable. 

It should be noted that when fixing the explicitly conformally flat form of the 
metrics (10), we thereby impose further restriction on the pseudogroup (5, namely, 
that the local diffeomorphisms be conformal maps. Fortunately, the classification 
problem of general two-dimensional metrics with respect to diffeomorphisms and the 
classification problem of metrics (10) with respect to conformal maps are identi- 
cal. Using the Wcyl-Lewis-Papapetrou coordinates r,s,w,f substantially simplifies 
computations on J°°Tr. 



3.5 A solution to the equivalence problem 

In this section we show how to reconstuct a generic metric g up to isometry from 
the knowledge of functional dependences between certain eight invariants computed 
for g. This provides a solution to the local equivalence problem as a particular case 
of the "principle of n invariants" [1, Ch. 7, § 4.3]. Let p, q be two independent scalar 
invariants, which p, q remain functionally independent when evaluated on a pseudo- 
Riemannian space g (genericity) . Then the metric g is locally uniquely determined 
by the values of six components gij{p, q) with respect to coordinates p, q; these com- 
ponents are themselves scalar invariants. Now, two metrics g, g' are locally equivalent 
if and only if the corresponding invariants gij{p, q),g,ij{p\ l') have equal expressions 
in terms of invariants p,q and p',q', respectively. That expressions for Eij{p,q) can 
be large is not a problem, since it is p, q which must be simple enough for the pro- 
cedure to be computationally tractable. In practice one usually finds a simpler set 
of invariants than gij{p,q), which is then suitable for storing in a database. Let g 
be fixed in the sequel. In what follows we interpret invariants as objects on S (cf. 
Remark 1). 

Proposition 4 Let p, q be two functionally independent invariants from the list Cp, 
C^, Q^, Q-y (on S). Then knowing the values of the remaining two invariants and 
the invariants Xp,Yp,Xq,Yq in terms ofp,q is sufficient for recovering the metric 
g G 

Recall that Xp, Yp, Xq, Yq are just four of the numerous second-order invariants 
of g discussed in Section 3.3. 

Proof. Once we know the values Xp, Yp, Xq, Yq in terms of p, q, we know the 
invariant frame X, Y in terms of coordinates p, q: 

d d d d 

X = Xp—+Xq—, Y = Yp—+Yq—. 
op oq op oq 

Knowing also the remaining two values of Cp, C^, Q^, Qy in terms of coordinates p, q, 
we can apply formulas given at the end of Section 3.2 to compute g, p, x, and hence 
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also 7, in terms of coordinates p, q. In particular, g is fully and uniquely restored in 
this way. 

To restore the remaining three components hij = we start with the 

determinant x, which is determined from g{X, ■ ) = dx/x up to a constant multiplier. 
In terms of Weyl-Lewis-Papapetrou parameters r,w,s, we have x = ^r^, hence r is 
also restored up to a constant multiplier. To obtain w,s, we distinguish two cases. 

If Q-y / 0, then 7 is the Cosgrove metric of constant curvature, i.e., a space 
form. As is well known, space forms are uniquely determined by their dimension and 
signature (see. e.g., Eisenhart [16, Ch. I, § 10]). Therefore, 7 is isometric to the form 
meaning that w, s can be restored up to isometry of 7 (isometrics 
of 7 reflect the freedom of determination of w and s by Proposition 3). This last step 
has been elaborated by Cosgrove [12] and reduced to solution of Appel equations. 

If Qj = 0, then by (9) we have 



s,i w^i 

S,2 W,2 



0. 



Therefore w, s are functionally dependent, i.e., s = s{z),w = w{z) for some function 
z on S, and 



7 = ^ ^2"^ (dz) 
s 



where prime denotes differentiation with respect to z. Assuming s non-constant, 
upon identification z = s we have 7 = (1 ^ Wg){ds/s)'^, which determines w{s) 
up to an additive constant. Assuming s = const and w 7^ const, the identification 
z = w leads to 7 = ^{1/ s'^){dw)'^ , which determines w up to an additive constant 
again. Assuming s = const, w = const would lead to = = and = \Cp, 
contradicting the assumptions of the proposition. □ 



3.6 Relation to \]/2 

The first-order invariants Cp, Cy., Q^., Q-y are not necessarily invariant with respect to 
general diffeomorphisms, since by very construction they depend on the uniqueness 
property of the subalgebra ^. In this section we relate Q-y to Im ^'2 and to Re ^'2, 
where 'I'2 is one of the Petrov invariants. 

Consider the Weyl tensor Cabcd = R-iema^ed - ga[cR-iCd]t + g6[cR'icd]a + 5R-ga[cgd]& 
in terms of adapted coordinates ,t'^ . Its nonzero components are 

C1234 = 5 V Q-y det h det g , 

C1212 = |(i? + i Alnx + ^ Cp - C^) det 5, 

C3434 = I {R + I Alnx + ICp - C^) det h, 

components of the form Cikji = 1,2 and k,l = 3, 4), which are too involved to be 
presented here, and components obtained from these through the symmetries Cabcd = 
—Cbacd = —Cabdc = Ccdab- Due to its Symmetries, the Weyl tensor can be thought 
of as acting on bivectors. The corresponding eigenvalue problem lies in the heart of 
the Petrov classification the Karlhede classification is a refinement of. To determine 
the Petrov type of a Lorentz metric, one examines the multiplicities of roots of the 
quartic '^ox^-4^^ix^y+Q'^2X^y'^-4^^3xy^+'^4y^, where *o = Cabcdk"'m''k''m'^, *i = 
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are calculated with respect to an arbitrary complex null tetrad {k,l,m,m,). Routine 
computation shows that ^'i = and ^'3 = (hence metrics in T cannot be of Petrov 
type ///) and then the quartic in question becomes ^'o-'^^ + 6^2X^Y^ + '^^Y'^. A 
metric is algebraically special if and only if ^'o^'4(9^2 ~ ^0^4) = 0. Moreover, ^'o^'4 
and ^'2 are invariants with respect to general diffeomorphisms. Of special interest in 
this paper is the invariant ^2 and its real and imaginary part. By routine computation 

Re^2 = ^C^-lR+^,K ^^^^ 

Im*2 = i V-Q-y 

(recall that Q-y < in case of Lorentz metrics). Thus, Im^2 is of the first order. 



3.7 Genericity 

It is the aim of this section to review the genericity assumptions we made in the 

course of exposition and express them in terms of scalar invariants. Here is the full 
list of the assumptions (additionally to det g and x = det h 0, which follow 
from the nondegeneracy of g) : 

1. The algebra of Killing vectors has a unique two-dimensional commutative sub- 
algebra M (Section 2); 

2. Cp^O (Section 3.2); 

3. There are at least two functionally independent invariants among Cp, C^, Q^, 
Qj (Proposition 4). 

The first of these assumptions, unlike the others, is not formulated in terms of 
invariants. Moreover, it cannot be formulated in terms of invariants that themselves 
depend on the uniqueness property of A. Yet we can improve the situation with the 
help of a classical result valid without restriction to the class T. 

Proposition 5 Let g be a general metric (possibly outside class for which there 
exist two functionally independent scalar invariants with respect to the general diffeo- 
morphisms. Then the distribution generated by Killing vectors of g is of dimension 
at most two. 

Proof. More general results are proved in Kerr [25, 26]. A simple proof goes as 
follows. Assume that ^1, ^2, ai'c three Killing vectors of g such that the distribution 
S they generate is of dimension three. Let p, q be scalar differential invariants. Then 
dp vanishes on H, since = ^ip = dp{^i), and similarly dq vanishes on H. But M. has 
only one dimension more than that of S, whence dpAdq = 0, proving that invariants 
p, q are functionally dependent. □ 

Turning back to metrics of class T, our Q-y is invariant with respect to general 
diffeomorphisms in view of results of the preceding section. The same is true about 
in case of vacuum Einstein metrics. 

Now we are left with the question when there exists a third Killing vector lying 
in the distribution H generated by ^. 

Proposition 6 The algebra of Killing vectors spanning the distribution H is three- 
dimensional if and only if coefficients hij are constant multiples of each other and 
two-dimensional otherwise. 
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Proof. In the adapted eoordinates we have = d/dv} , ^2 = djdv? ^ ^3 = d/du^ + 
d/du^ for some functions /* of t^,t^,u^,u^. The Kilhng equations ^3g = reduce 

to 

hik^ = 0, i,j = ^,2, (12) 



hju — T + huj — r = 0, 7 = 1,2. (13) 

Since det/i / 0, equations (12) imply that /* do not depend on t^,t^. 

Functions f^{u^, u^), P{u^,u^) now satisfy the remaining equations (13) and also 
derivatives of equations (13) with respect to t^,t^. The resulting nine homogeneous 
linear equations can be viewed as an algebraic system in four unknowns dP /dvf' , 
i,k = 1, 2. Its properties depend on the values of six 4x4 minors Dj, j = 1, . . .,6, 
formed by the rows corresponding to the three equations (13) plus one of the six 
derivatives thereof. These minors are 

Di = 8(/iii/ii2,i — /ii2^ii,i) det h, D4 = 8(/iii/ii2,2 — /ii2/iii,2) det h, 

D2 = 4(/lii/l22,l - /i22^11,l) det h, D5 = 4(/lii/l22,2 - /i22^11,2) dct h, 

D3 = 8(/ii2^22,i - /i22^i2,i) det h, Dq = 8(/ii2/i22,2 " /i22^i2,2) det h. 

If at least one of them is nonzero, then the algebraic system (13) has only the zero 
solution, hence are constants and ^3 belongs to ^. 

If, conversely, = for all j = 1, . . ., 6, then 

hij = Hcij (14) 

where H{t^,t^) is a nonzero scalar function and Cy is a nonsingular symmetric con- 
stant matrix (meaning that coefficients hij are constant multiples of each other). 
Under (14), system (13) has a unique solution 

d d 
ou ou 

modulo A and up to a constant multiple. □ 

By Proposition 6 we have either ^3 G .ft or (14). In the latter case = ACp, 
= 0, Qj = 0, contradicting assumption 3. 



3.8 Vacuum case 

We consider Einstein's vacuum equations with cosmological constant 

Ricij - I Rgij + Agij = (15) 

where A denotes the cosmological constant. 

Proposition 7 Let 7^ 4(5^. Let p = Cp, let q be an arbitrary invariant from the 
listC^, Q^, Qry, such thatp,q are functionally independent (onS). Then knowing the 
values of the remaining two invariants in terms of p, q are sufficient for recovering 
a generic metric g E % in the vacuum case. 
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Proof. Using the Weyl-Lewis-Papapetrou parametrization, we consider the system 
consisting of six Einstein's vacuum equations with cosmological constant (15) and 
eight equations for XCp, XC-^, XQ^., XQ^, YCp, YC^-, YQ^, YQ^. This system of 14 
equations is hnear in second derivatives of the four functions f,s,r,w. The number 
of second derivatives is ten {fu and fzz are missing). By applying linear algebra 
we obtain four nontrivial linear dependences among the second order invariants and 
components of Einstein's vacuum equations with cosmological constant. Rewritten 
in invariant form, they are 



YCp = -8/1, 

XCp = 8(Q^ - Q^) + Cp{C^ -ICp-AL), 

^ CjIiXC^ — ^ I\XQ^ + ^ l2YQ^ — ^ I^YC^ — I4YQ^ 

= h[h{C^ - i + 4A) - ^ C^Cpil Cp + 4A)], 
\ eCjIiYC^ — ^ eliYQ^ — ^ I2XQ^ + | I^XC-^ + I^XQ^^ 

= h{{\ Cp + 4A)(Q^ - Q^) + 2/3) - I hCpil Cp + 4A) 



(minus in elliptic case and plus in hyperbolic case). 

The derivation of these relations is tedious but straightforward and relies on 
inversion of a matrix with determinant (7^/4. While Cp 7^ is nonzero by assumptions 
(Section 3.7), /4 7^ is an additional gcncricity assumption. 

To finish the proof, we express Xp, Yp, where p = Cp, from the first and the 
second equation and Xq,Yq, where q is one of C^, Q^, Q^, from the third and the 
fourth equation. By Proposition 4 this is sufficient for recovering the metric g G T. 



For instance, if g = Q^, then expressing Xq,Yq is possible if and only if I4 7^ 0. 
Proposition 8 In the vacuum case, Re ^'2 is of the first order. 

Proof. On one hand the trace of (Ric|5)ij is 



where we denoted 




e = ±1 



□ 
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r^^'(Ric|5)ij =R-2A 



(16) 



by (15). On the other hand the trace of (5) is 



(17) 



Prom (16), (17) and by (6) we have 
a = - i Cp - 4A. 



(18) 
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Thus, Cv is a first order invariant in the vacuum case. The trace of (15) is 
R = 4A. 

Therefore, in the vacuum case we have 

R = - \C^. 
by (16), (17), (18) and (19). Then 

Re^-s = \C^ + \^ 
by (11). □ 

The other Petrov invariant ^^0^4 is of second order even in vacuum. 



(19) 
(20) 



4 Example 

Example 1 We ihustrate our results by classifying the Kerr-NUT-(anti)de Sitter 
space-time, attributed to Plebahski and Dcmiahski [37, 39], who derived it a special 

case of a seven-parameter family of Petrov type D solutions with nonzero cosmologi- 
cal constant A. Denoting the nonignorable coordinates as ti,t2, the metric coefficients 
are 



/ill = 



Q 



tl + tl 



ti+tr 



hi2 = 



Ptl 
tl + tl 



+ 



Qtl 

tl + tV 



h-22 = 



tl + tl 



Qtj 

tl + tV 



and 



where 



511 



tl + tl 



912 = 0, 



522 



tl + tl 



Q 



P={A^- tl){l + I Atl) + 2Lti, Q = {A'^ + tl){l - I Atl) - 2Mt2, 

M, L, A, A being parameters of the metric. While M is the mass of the source, the 
physical meaning of L,A is not immediately obvious. Griffiths and Podolsky [18] 
gave the following identification 

^,^^^,_^,^(a^ + 3aA + 3 



{a^-f)A + 3 ' 



L = I 



3/2)A + 3 



(a2-/2)A + 3 J' 



where a is the angTilar momentum and I is the NUT parameter. 

The four basic first-order invariants Cp,Cy-,Q^,Q^, can be routinely computed. 
Two of them turn out to be particularly simple, hence we choose them as the inde- 
pendent invariants p, q: 



P = Q'y = - 
q = C^=4 



Lt2{tl - 3tl) - MtiiStl - tl) 



HI + tl)' 

Mt2{tl - 2,tl) + LtiiStl - tl) 

2\3 



(tf + tl) 



-1a. 

3 



(21) 
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Components of the invariant frame are easily seen to be 
2 2Atl - (A^2 - 3)ii - 3L 



tl + 4 



^2 _ _ 2 2A4 + {AA^ - 3)t2 + 3M 



tl + 4 



Y^=sgnQ 



= -sgnP 



X' 



X\ 



Expressions for the other invariants needed for recognizing equivalence are quite 
large. They are schematically given in the following table: 



Invariant 


Numerator 


Denominator 




Pl0,6,6 


{tl + tl)PQ 


Qx 


Pl7,14,14 


{tl + tlfPQ 


XQy 


P3,3,3P7,7,7 


{tl + tlf 


XC^ 


P7,7,7 


(tl + tlf 


YQ^ 


P3,3,3 Pll,8,8 






Pll,8,8 


(t? + ti)5pV2Ql/2 



Here Pn,ni,n2 is a substitute for a polynomial p(ti,t2) of total degree n, degree ni in 
the indeterminate ti and degree n2 in the indeterminate t2 (e.g., the substitute for 
P would be P4,4,o)- 

Actual values of these invariants are suitable for storing in a database of exact 
solutions. Prior to that one can wish to express all of them in terms of p, q. This 
usually constitutes the main technical difficulty connected with equivalence analysis. 
In our example a help comes from the simple relation 



iiA+c^r 
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Inserted back to the denominators in (21), this allows us to find the following cubic 
equations for ti and t2 in terms of invariants p, q and parameters A, M, L: 



where 



4t? 

4tl 



I = 



3Iti + /4 
3/^2 - I- 



0, 
0, 
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{n+±Ay-iQ, 



1/3 



/± = 4 



M{C^ + jA)±2L^/^^ 



One easily checks that /4 7^ so that by Proposition 7 the metric is characterizable 
by expressions for C^, in terms of Cp, Q^. 

Example 2 In the special case of A = (Kerr NUT solution) we have L = I, 
P = —tl + + 2 Lti, Q = t\^ — 2 Mt2 quadratic in ti,t2, respectively, and the 
invariants simplify considerably: 
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In\'ariaiil 


Numeral or 


Dcnoniiualor 




P2,2,2 


{tl + tl)PQ 


Qx 


Pl0,8,8 


{tl + tlf PQ 




^•3,3,3^5,5,5 


{tl + tlf 


XC^ 


P5,5,5 


{tl + tlf 


YQ^ 


^3,3,3^*7,6,6 


{tl+tlf p^'^Qy' 


YC^ 


P7,6,6 


{tl + tlf py^Q^-" 



The same Kerr-NUT solution has been classified by Aman and Karlhede in [2]. 

These authors start with the metric [27, II. A], given in the Newman-Penrose tetrad 
formalism. Since this metric is not in the form (1), the results cannot be compared 
immediately. In the course of the Karlhede algorithm applied to this metric, 14 
curvature invariants are computed. These are complex quantities, whose real and 
imaginary parts are to be considered separately. This example confirms that a tai- 
lored procedure can lead to a significant reduction of the total number of computed 
invariants. 

On the other hand, all the Aman-Karlhede invariants already have (nearly) ex- 
plicit expression in terms of two complex conjugated quantities p,p*-, which them- 
selves are easily derivable from the unique complex curvature invariant of order zero. 
Thus the "last step" becomes rather trivial in this case. However, this may be at- 
tached to the particular way the Kerr-NUT metric (which is of Petrov type D) has 
been derived. 

5 Conclusions 

We have demonstrated that four first-order scalar differential invariants and two first- 
order invariant vector fields suffice to solve the local equivalence problem within the 
class T of pseudo-Riemannian metrics with two commuting orthogonally transitive 
vector fields. Results are restricted to generic metrics in T as identified in Sect. 3.7, 
leaving the special subclasses unexplored. The main open problems include: the de- 
scription of the special (nongeneric) subclasses exempt from consideration in this 
paper; the classification problem within these classes; the comparison to the Karl- 
hede algorithm; finding solutions to the Einstein equations in terms of first-order 
invariants (cf. [7] and references therein). 
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